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Summary 

Soil erosion by water, wind and tillage affects both agriculture and the 
natural environment. Interest in fallout radionuclides such as Sr and 
137Cs to study medium term soil erosion (40 yrs) increased in the early 
1990’s. The mass balance approach used to estimate long-term erosion 
and sedimentation rates is presented. A modification to improve the 
model for areas which are subject to deposition of eroded soil is 
presented. As an example, data from a site close to Berlin are presented. 
Using our improved mass balance model, !*’Cs derived estimations of 
erosion and deposition rates are presented. 


1. Introduction 


Studying long-term erosion and sedimentation processes by using man- 
made and natural radioisotopes is an innovative approach in soil science, 
which has developed over the past 30 years. Fallout '%’Cs is the 
radionuclide most extensively been used to provide independent 
measurements of soil-erosion and sediment-deposition rates and patterns 
[1, 2, 3, 4]. 


1.1 Erosion measurements using ‘*’Cs 


Caesium-137 from atmospheric nuclear-weapons tests in the 1950s and 
1960s (Fig.1) is a unique tracer of erosion and sedimentation, since there 
are no natural sources of '%’Cs. Events such as the Chernobyl accident in 
April 1986 caused regional dispersal of '%’Cs that may affect the 
deposition budget at a specific site. 
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Figure 1: Mean annual !°’Cs fallout deposition rates in Germany (50° N- 
60° N) (modified from [8]) 
137Cs can be easily measured by gamma spectroscopy. Using !°’Cs is a fast 
and cheap method to study erosion-deposition processes compared to 
traditional methods like silt bags. 


Because of its long half life, '*’Cs may migrate into the soil. Its 
distribution in soils is affected by spatial variations in deposition rates 
and by land use. An undisturbed soil will exhibit higher radionuclide 
activities near the soil surface, which reflects their surficial input and 
slow downward transport. Plowing homogenizes '°’Cs within the plowed 
layer (Figure 2). The concentration of fallout '°’Cs in surface horizons and 
its slow transport into the soil have allowed using this radioisotope to 
estimate rates of soil erosion [4, 5, 6]. 
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Figure 2: Generalized sketch illustrating the distributions of '%’Cs in soils 
under tillage and undisturbed conditions. 


2. Erosion study at Young Moraine regions of Germany 


Soil redistribution rates were estimated at a real site in the Young 
Moraine region of North-East Germany. The sampling was performed on 
16” April 2008 at Basedow, a cultivated area north of Berlin (53° 34’ N & 
13° 80’ E). Figure 3a shows the transect used for sampling. Sites which 
were not affected by soil loss or deposition and thus reflect atmospheric 
inputs of '°’Cs and its decay were taken as reference sites [7]. The 
sectioned soil cores were collected from depths of 0-10 cm in increments 
of 2 cm and from depths below 10 cm in increments of 5 cm. Samples 
from 15 measurement sites and 2 reference sites were taken. The soil 
was dried at 105°C, sieved, filled in aluminium bottles and analysed by 
gamma spectroscopy. Figure 3 shows the measured '*’Cs inventories 
along the transect under study. 


2.1 Conventional mass balance erosion/deposition model 


The standard approach to erosion assessment using '°’Cs is described and 
discussed in detail elsewhere [1, 2, 3, 4]. The inventory of '°’Cs at the 
measurement points (A:) is compared with a reference inventory (Aref) 
taken at areas not subject to soil erosion or deposition. If A: < Ares erosion 
is assumed to have occurred; A: > Aree is assumed to result from soil 
deposition. The '°’Cs concentration may vary due to cultivation practices 
such as ploughing. The rate of soil erosion and deposition can be 
estimated as: 
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Where Yer is the erosion rate (t hat yr‘), X the reduction in total '°’Cs 
Aw A 


inventory ( ), d the ploughing depth (m), ta is the time step ( 


Ar 
ta =b +n-At,n=0 set to t =1963) and p, the soil’s bulk density (kg m*). 
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Figure 3: Top: study transect at Basedow; bottom: Activity distribution of 
137Cs along the transect; the dotted band gives the 95% uncertainty 
interval of the reference activity, error bars denote sample standard 
deviations. 


The eroded soil is carried by the force of water and wind and is 
deposited. The mean of the concentration (Bq m°) of '*’Cs deposited 
annually from the upslope areas, C,(¢), is given by eqn. (2), the 
deposition rate at year n is calculated by eqn. (3): 
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Here C;,(4)denotes the '°’Cs concentration mobilized from eroded area S; 


(Bg m3), AZ the height of soil eroded during a time step tn, Si the slope 


area represented by point i (m°), P is the particle size correction factor, K 
is the number of eroding areas and N is the year of sampling. 


2.2 Improved Mass Balance Deposition Model 


The basic assumption of the mass balance approach is that the mass of 
soil eroded must be deposited in the area studied. The deposition sites 
thus reflect the amount of dislocated soil. The initial concentration of soil 
containing '°’Cs is Co. In the following years a constant layer, AZa, of soil 
eroded from upslope areas is deposited in the area down the slope. As a 
consequence the depth, to which the soil (and the '°’Cs sorbed to it) is 
homogenized due to ploughing, changes annually. This effect that is 
neglected in eqn. (3) is taken into account in the following. 


The deposited concentration Cit) (Bg m°) is then estimated as 
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where Ca (Bq m°) is the '*’Cs concentration of the eroded soil given by 
eqn. (2). 


The deposition rate (t ha! yr’) is estimated as 
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Figure 4: '°’Cs concentration at the depositional site predicted for 2008 
by eqn. (4) assuming different erosion rates. 


As a first step to explore the effect of the improvement supposed in our 
soil deposition model given by Eqn. (4), numerical simulations were 
performed. Figure 4 shows depth distributions of deposited '°’Cs at the 


depositional site for four erosion rates as predicted for the year 2008 
using the deposition model of Eqn. (4). In Fig.4 the maximum of the '°’Cs 
concentration corresponds to the '°’Cs fallout in 1963, and thereafter its 
concentration decreases each year, since at the eroding site each 
ploughing leads to mixing with some soil without Cs. The homogeneous 
distribution of '°’Cs in the upper 30 cm results from ploughing. 


3. Estimation of deposition rates using the improved mass 
balance model 


The !8’Cs mass balance model of Eqn. (1) was used for calculating erosion 
rates at the two segments of our experimental site shown in Fig. 3. The 
deposition rates calculated by equations (3) and (5) are given in Table 1. 
The differences between deposition and erosion rates (Apg) show that the 
deposition rates calculated by eqn. (3) at both sites violate the mass 
balance of eroded and deposited soil. 
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Figure 5: Simulations of deposited '°’Cs concentrations at the two 
segments (top: site 1, bottom: site 2 of Fig. 3) based on egn. (3) (* signs) 
and (5) (+ signs), respectively. 


Fig. 5 show that eqn. (3) overestimates deposited concentrations and, as 
a consequence, underestimates deposition rates compared to the 
improved model of egn. (5). This result from including the effect of annual 
reduction of '°’Cs by the fraction no longer ploughed due to the increase 
in soil deposited. 


Table1: Soil budget for Basedow site derived from two mass balance 
deposition models; Ap: gives the difference between deposition and 
erosion rates 


Erosion Deposition rates calculated (t ha" yr' ) 
rates Egqn.(3) ADE Eqn.(5) ADE 
(t ha yr! 
) 
Site 1 -18 7 11 17 1 
Site 2 -44 26 18 46 2 
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